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The understanding of the interaction of nucleons in nuclear and neutron-rich matter at non- 
zero temperature is important for a variety of applications ranging from heavy-ion collisions to 
nuclear astrophysics. In this paper we apply the Dirac-Brueckner-Hartree-Fock method along with 
the Bonn B potential to predict single-particle properties in symmetric nuclear matter and pure 
neutron matter at finite temperature. It is found that temperature effects are generally small but 
can be significant at low density and momentum. 



I. INTRODUCTION 

The nuclear equation of state (EoS) at finite temper- 
ature is of fundamental importance for heavy-ion (HI) 
physics [J as well as nuclear astrophysics, particu- 
larly in the final stages in the evolution of a supernova. 
Knowledge of the finite-temperature EoS can be of great 
help in the interpretation of experiments aimed at iden- 
tifying a liquid-gas phase transition. Presently, find- 
ings concerning such transition are very model dependent 

When other aspects, such as spin asymmetries of nu- 
clear/neutron matter are included as well, conclusions 
are even more contradictory. For instance, the influence 
of finite temperature on the manifestation of ferromag- 
netic instabilities is an unsettled question. In Ref. 10], 
phenomenological Skyrme-type interactions are used in a 
Hartree-Fock scheme. It is found that the critical density 
for ferromagnetism decreases with temperature. On the 
other hand, in Ref. [11] the authors report no indication 
of ferromagnetic instability at any density or temperature 
based on the Brueckner-Hartree-Fock (BHF) approxima- 
tion and the Argonne V18 nucleon-nucleon (NN) interac- 
tion ^12j . The properties of spin-polarized neutron mat- 
ter at finite temperature are studied in Ref. [l^ with 
two different parameterizations of the Gogny interaction. 
The results show two qualitatively different behaviors for 
the two parameterizations. The reasons for these discrep- 
ancies must be carefully studied and their origin under- 
stood in terms of specific features of the nuclear force 
and/or the chosen many-body framework. 

Previous work on the temperature dependence of the 
EoS includes the calculations by Baldo and Ferreira who 
used the Bloch-De Dominicis diagrammatic expansion 
, Brueckner-Hartree-Fock calculations with and with- 
outh 3BF 0, and the predictions of Ref. [IB] based on 
the Green's function method. The entropy per particle in 
symmetric nuclear matter has been studied in Ref. [l3] 
within the self-consistent Green's function (SCGF) ap- 
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proach, where both particle-particle and hole-hole scat- 
terings are included. The SCGF framework allows direct 
access to the single-particle spectral function and thus to 
all the one-body properties of the system. A most re- 
cent work by Rios et al. [iTj addresses hot neutron mat- 
ter within the same approch and performs a comparison 
with other models. Hot asymmetric matter and /3-stable 
matter have been studied by Moustakidis et al. (Tsl [l9j 
using temperature and momentum dependent effective 
interactions. 

The study of the many EoS-related aspects in both 
symmetric and neutron matter starting from realistic NN 
forces and within a microscopic model is still a consider- 
able challenge. It is the purpose of this paper to report 
the first part of a comprehensive study of temperature 
dependence of nuclear and neutron-rich matter proper- 
ties based on the Dirac-Brueckner-Hartree-Fock (DBHF) 
approach. An earlier calculation with the DBHF method 
can be found in Ref. However, our work will go 

beyond those predictions, as we plan to address addi- 
tional aspects such as (spin/isospin) asymmetric matter 
and temperature dependence of in-medium effective cross 
sections. Previously, we have confronted isospin asymme- 
tries [10, [2l|, [l^] as well as spin asymmetries [1^ effects 
on the equation of state of cold matter. The inclusion 
of temperature dependence and the simultaneous con- 
sideration of all of the above mechanisms will make our 
microscopic predictions more broadly useful and capable 
to reach out to the properties of the hot environment 
present in the latest stage of a supernova collapse or in 
the collision of heavy nuclei at intermediate energies. 

Next, after a review of the main aspects of the formal- 
ism, we will show and discuss predictions of the chem- 
ical potential and single-particle properties in symmet- 
ric and pure neutron matter. We will concentrate here 
on the properties of the single-particle within the nu- 
clear medium for the following reason: although the nu- 
clear/neutron matter energy density is certainly an im- 
portant quantity, the single-particle interaction and its 
temperature dependence, which determine the one-body 
properties in the medium, are perhaps more relevant for 
non-equilibrium processes such as relativistic heavy-ion 
collisions. 
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Hot nuclear and neutron matter are infinite fermionic 
systems where only the strong interactions among nucle- 
ons are taken into account. The temperatures typically 
considered are of a few tens of MeV, relatively small on 
the scale of nuclear energies. For instance, at densities 
near the saturation density of nuclear matter, the free 
Fermi energy, e^, is approximately 40-50 MeV, and thus 
a temperature of 20 MeV is still somewhat low, with 
T/ep 0.5. 



II. FORMALISM 

Within the DBHF method, the interactions of the nu- 
cleons with the nuclear medium are expressed as self- 
energy corrections to the nucleon propagator. That is, 
the nucleons are regarded as "dressed" particles, essen- 
tially a gas of non-interacting quasi-fermions. The be- 
havior of the dressed nucleon is determined by the effec- 
tive nucleon propagator, which obeys the Dyson equa- 
tion. Relativistic effects lead to an intrinsically density- 
dependent interaction which is consistent with the con- 
tribution from three-body forces (TBF) typically em- 
ployed in non-relativistic approaches. The advantage of 
the DBHF approximation is the absence of phenomeno- 
logical TBF to be extrapolated at higher densities from 
their values determined through observables at normal 
density. 

In the quasi-particle approximation, the transition to 
the temperature-dependent case is introduced by replac- 
ing the zero-temperature occupation number with its 
finite-temperature counterpart, namely the Fermi-Dirac 
occupation density. In fact, it can be shown that the 
zero-temperature terms in the Bethe-Goldstone expan- 
sion where temperature is introduced in the occupation 
number only are the dominant ones thus justifying 
this simplified procedure. More specifically, one replaces 



n{k,p) = 



1 \ik<kF 
otherwise, 



with the Fermi-Dirac distribution 
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Here T is the temperature in MeV, e{k,p,T) the single- 
particle energy, and ^ the chemical potential to be deter- 
mined. The angle-averaged Pauli operator is evaluated 
numerically. It's given by 



1 

Q{q, P,p,T)^- J d{cos e){l - npD 



(fci,p,r))x 



X (1 - nFD{k2,p,T)) 



(3) 



for two nucleons having momenta ki and fc2 , with relative 

and total momentum given hy q — and P — fci-f fc2, 

respectively. 



The single-particle energy, e{k,p,T), is now tempera- 
ture dependent. It can be obtained self-consistently with 
the Dirac states following the same procedure as used in 
the zero temperature case, but including Eq. (2) in the 
calculation of the single-particle potential. At each iter- 
ation of the self-consistent calculation, the normalization 
condition 



P 



1 f°° 
^TT^ / nFD{k,p,T)d'k 



(4) 



allows to extract the microscopic chemical potential, 
p{p,T). D is a degeneracy factor, equal to 4 for sym- 
metric unpolarized nuclear matter or 2 for unpolarized 
neutron matter. 

In close analogy with the T = case, the single- 
particle potential and the self-consistent nucleon-nucleon 
G-matrix are related by 



U{k,p,T)^ 



(2J + 1)(2J+1) 
(2t+ l)(2s + 1) 



X/ nFD{k',p,T)Gii^''''{q{k,k'),P{k,k'))d^k', (5) 
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where /, L, S, and J are the NN system quantum num- 
bers while s and t are the single-particle spin and isospin, 
respectively. Integrations over the Fermi sea are per- 
formed using a cutoff adjusted to the density. We found 
a value of 6 fm to be sufficient for the highest densities 
and temperatures considered here. 

The single-particle energy is the sum of potential ([/) 
and kinetic {KE) energy contributions. 



e{k,p,T)^Uik,p,T)+KE{k,p,T). 



(6) 



It is parametrized using the same ansatz as in the T — 
case (24i] for a spin-symmetric, rotationally invariant 
system: 



e(fc,p,T) - ^k^ + {m* f + Uv, 



(7) 



where C/y is the time-like part of the vector potential 
and m* — m -\- Us, with Us the scalar potential, now 
both density and teperature dependent. This quantity 
plays an important role in temperature and momentum- 
dependent transport models of heavy-ion collisions. 

Once a self-consistent solution is obtained for the 
single-particle potential, and, simultaneously, for the 
chemical potential, the entropy/particle can be calcu- 
lated, along with all thermodynamic quantities. In the 
mean- field approximation, the entropy /particle is given 

by 



S 



D 



(27r)3 



[npoik, p,T) In npoik, p,T)- 



+ (1 - npoik, p, T)) In (1 - npoik, p, T))] £k, (8) 

that is, it has the same functional form as in a nonintcr- 
acting system. 
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FIG. 1: (color online) The chemical potential in symmetric 
matter as a function of the Fermi momentum at various tem- 
peratures from T = to T = 30 MeV in steps of 5 MeV. The 
chemical potential decreases with temperature. 
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FIG. 4: (color online) The single-particle energy minus the 
rest mass in symmetric matter parametrized as in Eq. (7) vs. 
the momentum at different temperatures: T = (solid black) ; 
T = 10 (dashed green); T = 20 (dotted red). The left and 
right panels correspond to Fermi momenta equal to 0.9 fm~^ 
and 1.3 fm~^, respectively. 



III. RESULTS 




FIG. 2: (color online) The single-particle potential in sym- 
metric matter as a function of the momentum for three dif- 
ferent temperatures: T — (solid black); T = 10 (dashed 
green); T = 20 (dotted red). The left and right panels corre- 
spond to Fermi momenta equal to 0.9 fm~^ and 1.3 fm~^, 
respectively. 
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FIG. 3: (color online) Effective mass in symmetric matter as 
a function of the Fermi momentum and for different temper- 
atures: T = (solid black); T = 10 (dashed blue); T = 15 
(dash-dotted green); T = 20 (dotted red). 



A. Symmetric Nuclear Matter 

We begin by showing the pre(iicted chemical potential 
in symmetric nuclear matter (SNM), see Fig. 1. Each 
curve is an isotherm starting from zero temperature and 
going up to T = 30 MeV in steps of 5 MeV (/i goes down 
with increasing temperature). The effect of temperature 
is definitely larger at the lower densities and increases in 
size with increasing temperature. Our predictions are in 
fair qualitative agreement with those shown in Ref. p^ . 
with or without the contribution of TBF, which do not 
seem to have a major effect on the chemical potential. 

Next, we show the single-particle potential in SNM, 
U{k,p,T), see Fig. 2. The effect of temperature is much 
more pronounced at low density (compare left and right 
panels) and low momenta. The general tendency is to 
turn slightly more repulsive with increasing temperature, 
although this trend becomes clear only at the higher 
temperature. The increased repulsion is the result of a 
combination of effects. Temperature "smears out" the 
step function distribution, Eq. (1), so that the interac- 
tion probability increases (decreases) at low (high) mo- 
menta due to the smaller (larger) occupation probability 
as compared to the T = case. Although Pauli blocking 
is generally reduced by the onset of temperature (which 
would suggest increased attraction among the particles), 
the integral in Eq. (5) receives contributions from G- 
matrix elements at higher momenta (as compared to the 
zero temperature case), and such contributions tend to 
be repulsive. In the end, we observe a net effect that 
is repulsive, except at the lowest temperatures. Consis- 
tently with Fig. 2, the temperature dependence of the 
effective mass is generally small, see Fig. 3, and more no- 
ticeable at low density. Thus, in the range of densities 
and temperatures considered here, we expect only minor 
effects on the in-medium cross sections, whose behavior 
is essentially dominated by the effective mass. A slightly 
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FIG. 5: (color online) The entropy /particle in SNM as a func- 
tion of density and for increasing temperatures of T = 5, 10, 
15, and 20 MeV. 
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FIG. 6: The entropy /particle in SNM as a function of tem- 
perature at a density corresponding to a Fermi momentum of 
1.3 fm-^. 



repulsive temperature effect on the DBHF effective mass 
was found in Ref. |8i] as well. 

The single- particle energy, not including the rest mass, 
is shown in Fig. 4 at fixed density and for different tem- 
peratures. Fig. 4 confirms that the single-particle interac- 
tion is noticeably impacted only at the lowest momenta. 
In Ref. Q the equivalent Schroedinger optical potentail 
was constructed from the Dirac self-energies and found 
to be remarkably insensitive to temperature. 

We conclude this set of results by showing a macro- 
scopic quantity, namely the entropy, which is a measure 
of thermal disorder. Entropy production in multifrag- 
mentation events in heavy-ion collisions is a crucial quan- 
tity in the determimation of the mass fragment distribu- 
tion. The entropy per particle is shown in Fig. 5 as a 
function of density and for various temperatures, and in 
Fig. 6 it is displayed as a function of temperature at a 
fixed density (close to saturation density). The entropy 
increases with temperature, as physically reasonable, and 
decreases substantially with density. At low T, it is ex- 
pected to approach a linear dependence due to the fact 
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FIG. 7: (color online) The chemical potential in neutron mat- 
ter as a function of the Fermi momentum at T = 0, 10, and 
20 MeV. The chemical potential decreases with temperature. 
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FIG. 8: (color online) The entropy /particle in NM as a func- 
tion of density at T = 10 (lower curve) and 20 MeV (upper 
curve) . 

that, for a Fermi liquid, the relation between S and T 
should be approximately 

S « ^NiT = 0)r, (9) 

in terms of the density of states at the Fermi surface. 

Although weak model dependence is not a general fea- 
ture of thermodynamic quantities, the authors of Ref. 
demonstrate that different approximations to the en- 
tropy, including the quasi-particle approximation in the 
temperature dependent BHF scheme, differ from each 
other by 10 to 20% at most. This may be due to cancela- 
tions in the difference between the single-particle energy 
and the chemical potential p^ . 

B. Neutron Matter 

We have done a similar study of hot neutron matter 
(NM) as well. The chemical potential in NM is shown 
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FIG. 9: (color online) The single-particle potential in neutron 
matter as a function of the momentum at different temper- 
atures: T = (solid black); T = 10 (dashed blue); T = 20 
(dash-dotted green); T — 30 (dotted red). The neutron mat- 
ter Fermi momentum is equal to 0.9 fm~^. 



in Fig. 7 as a function of the Fermi momentum. (Of 
course, tfie same Fermi momentum corresponds to half 
the density as compared to SNM.) The microscopic chem- 
ical potential, as obtained from the normalization condi- 
tion Eq. (4), was shown to be in good agreement with 
the macroscopic one, obtained from the bulk properties 
through the derivative of the free energy density IT]. 

The entropy in NM is displayed in Fig. 8. As for the 
case of SNM, discrepancies between predictions from dif- 
ferent models have been found to be small, especially 
those arising from the use of different NN potential mod- 
els In fact, our predictions (based on the Bonn 
B potential) are close to those shown in Ref. [l^ either 
with CDBONN ll] or Argonne V18 Concerning dif- 
ferent many-body approaches, there are indications that 
contributions to the entropy from dynamical correlations 
(which would fragment the quasiparticle peak) are small. 
Hence, different predictions tend to approach the "dy- 
namical quasiparticle" result [l7| . 

The temperature dependence of the single-particle in- 
teraction in hot NM is comparable to, or smaller than, 
the one encountered in the SNM case. We show a repre- 
sentative case in Fig. 9. As discussed in Section IIIA, we 
notice, from Fig. 9 and from the left panel in Fig. 2, that 
the potential tends to become deeper at first and then 
more shallow as the temperature increases. Also, tem- 
perature appears to "wash out" some of the structure in 
the potential. 



sion of Pauli blocking, which is most important at low 
momentum, some temperature-induced modification of 
(low-energy) in-medium cross sections could be expected. 
This is an interesting point to be explored further. We 
also discussed the entropy/particle. When compared 
with other predictions in the literature, our results con- 
firm the weak model dependence of this quantity. 

As usual, we adopt the microscopic approach for our 
nuclear matter calculations. Concerning our many-body 
method, we find DBHF to be a good starting point to 
look beyond the ground state of nuclear matter, which it 
describes successfully. The main strength of this method 
is its inherent ability to effectively incorporate crucial 
TBF contributions [24I yet avoiding the possibility of 
inconsistency between the parameters of the two- and 
three-body systems. 

We have focussed on the one-body properties of the 
system, which are of great relevance for the study of en- 
ergetic heavy-ion collision dynamics. A continuation of 
this study, including a close look at all thermodynamic 
quantities, will be the topic of a forthcoming work. 
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IV. SUMMARY AND CONCLUSIONS 



We used the Dirac-Brueckner-Hartree-Fock method 
extended to finite temperatures to predict single-particle 
properties in hot SNM and NM. For temperatures up 
to a few tens of MeV the effect of temperature is small 
except at low densities and momenta. Due to suppres- 
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